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Abstract 

We present the full thermodynamics of a fluid confined by an arbitrary external potential based on 
the virial expansion of the grand potential. The fluid may be classical or quantum and it is assumed 
that interatomic interactions are pairwise additive. We indicate how the appropriate "generalized" 
volume and pressure variables, that replace the usual volume and hydrostatic pressure, emerge for 
a given confining potential in the thermodynamic limit. A discussion of the physical meaning and 
of the measurement of these variables is presented. We emphasize that this treatment yields the 
correct equation of state of the fluid and we give its virial expansion. We propose an experiment 
to measure the heat capacity, so that with this quantity and the equation of state, the complete 
thermodynamics of the system may be extracted. As a corollary, we find that the so-called local 
density approximation for these systems follows in the thermodynamic limit, although we also point 
out that it cannot be used indiscriminately for all local variables. Along the text we discuss the 
relevance of these findings in the description of the currently confined ultracold gases. 
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largely stimulated first principles analysis of the physics of these inhomogeneous systems. In 
particular, the calculation of density profiles has turned out to be of fundamental importance 
since it is one of the main measurable quantities in the recent experiments. From the density 
profile several thermodynamic quantities can be obtained, such as the number of particles 
and the temperature. However, because the hydrostatic pressure is a local quantity and the 
volume of the system cannot be rigorously defined, the use of thermodynamics as a tool to 
analyze and characterize such trapped gases has been of limited application. For instance, 
important properties and quantities, such as the equation of state of the fluid or its heat 
capacities, are lacking. As we discuss here, those properties can actually be extracted from 
the knowledge of the density profile and other simple properties of the trapped gases. 

The main purpose of this article is to emphasize the fact that the thermodynamics of 
trapped systems must be reformulated in terms of the appropriate mechanical variables that, 
for lack of a better name, we shall call generalized pressure V and volume V. That is, the 
usual hydrostatic pressure p and the volume of a fluid contained in a vessel of rigid walls 
are no longer thermodynamic variables for a fluid confined by an external inhomogeneous 
field Vea;j(r). As mentioned above, in the presence of an arbitrary confining potential, the 
pressure of the fiuid becomes a local variable p = p{f) and the volume is strictly undefined. 
As we shall show there is a "new" and unique pair of variables V and V that replace the 



usual ones, p and V. We shall illustrate these variables with specific examples. A 
result has already been pointed out for a gas trapped in a harmonic potential 



though this 

15, y, m 



its use has been limited within the approximations that the gas is ideal or at best that 
the interactions may be treated a la Hartree-Fockjl^; nevertheless, incipient comparisons 



261] showing its potential 



have been made with experiments in ultracold Na gases [24 
usefulness. In any case, we believe the main point has not been fully appreciated, namely, 
the fact that a different set of thermodynamics variables must be used for a given external 
confining potential. We proceed here to fill this gap, deducing the virial expansion for 
the Grand Potential in the grand canonical ensemble for an arbitrary external potential. 



valid for classical and quantum fiuids, either ideal or with pairwise interatomic interactions. 
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Extensions to three- or higher-body interactions may be further considered. We discuss 
the physical meaning of the generahzed pressure and its relevance regarding the equation 
of state of the fluid, i.e. V = V{V/N,T) with N the number of molecules or atoms in the 
fluid and T the temperature. This finding should have immediate practical applications 
since all the thermodynamic properties, specially phase transitions, can be quantitatively 
described and certainly visualized in the corresponding phase diagram. We emphasize the 
fact that the generalized pressure can be very simply measured or calculated once the particle 
density profile p(r*) is known. This procedure has already been exploited in the analysis 
of experimental data of trapped Na ultracold gases in quadrupolar 2J] and in harmonic 



traps 
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261 ]. As a further step, we propose a simple and independent experiment that. 



in addition to the knowledge of the equation of state, allows for the determination of the 
heat capacity at constant generalized volume V, Cy. This experiment should be easily 
performed in the currently confined ultracold gases. To the best of our knowledge there are 
no measurements of such a heat capacity. It is a simple exercise to show that knowledge 
of the equation of state V = ViV /N,T) and the heat capacity Cv = C\;{N,V,T) suffices 
to know all the thermodynamics of a pure fiuid. We need not overemphasize the fact that 
the quantitative features of the latter quantities are direct consequences of the interatomic 
interactions and of the collective interactions of the fiuid. 

As a corollary of our analysis we shall show that the so called "local density approxima- 
tion" (LDA) follows within the appropriate thermodynamic limit of the confinin g po tentials. 
Our results are i n ag reement with rigorous proofs of LDA for classical systems 27|, 1281] and 
quantum systems 291] as well; these works have been largely overlooked in the current liter- 
ature of ultracold gases, but are very relevant since they show that, under the appropriate 
conditions, LDA is an exact procedure. The validity of LDA is expected, and certainly 
widely used, because in the trapped gases the main non-uniformity appears at macroscopic 
length scales due to the presence of the confining external field. That is, the latter must 
be "macroscopic" in order to trap a large number of particles. For macroscopic potentials 
this equivalently means that the energy-level separations of the external potential are much 
smaller than the typical atomic and collective excitations 

Although LDA has been used in essentially all current works on ultracold gases that yield 
density profiles, its validity is mostly assessed on these physical grounds. Thus, it is reassur- 
ing to find that the virial expansion when applied to confined quantum fiuids yields LDA as 
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well. As we shall discuss, however, care must be taken when using it. That is, LDA can be 
directly used for the grand potential density, the entropy density and the particle density, 
but not so for the internal energy and other free energy densities, among other possible 
thermodynamic quantities. That is, while it gives a procedure to calculate thermodynamic 
properties of a confined inhomogeneous fluid, it does not imply that the local states of the 
trapped fluid are thermodynamic states of the corresponding homogenous fluid. 

We proceed as follows. First, we describe the general system under study - an interacting 
fluid confined by an external, inhomogeneous potential Vea;((r) - and we provide several spe- 
cific examples to introduce the generalized volume V and the corresponding thermodynamic 
limit. By analyzing the system in the grand canonical ensemble, we derive the general virial 
expansion for the Grand Potential. We give the minimum conditions that must be obeyed by 
the interatomic interaction potentials for the expansion to be valid. We discuss the physical 
interpretation of the generalized pressure V and show how the equation of state may be 
found from the knowledge of the density profile; we propose a simple experiment to measure 
the heat capacity at constant V in the current experiments of ultracold gases. As mentioned, 
LDA follows as a corollary and we discuss some of its consequences and practical uses. 



II. GENERALIZED VOLUME AND THE THERMODYNAMIC LIMIT. 



The system consists of identical atoms or particles of mass m with Hamiltonian 



N -2 N 

Hn = T.^ + I: + E VeAri)- (1) 

i=l i<j i=l 

For simplicity of presentation we have assumed additive pairwise potentials but, as can be 
seen below, the analysis may be extended to arbitrary interatomic interactions. The external 
potential Vexti'f^ is responsible for confining the system. To serve this purpose, it should 
have at least one minimum and must obey that Vexti'f^) oo for |r| ^ oo. For rigid- wall 
containers it is costumary not to write down the potential. Here, we include it as Vext{r) = 
if r is within the volume V enclosed by the rigid walls and = cxd if r is outside of it. 

Typical examples of traps of atomic gases are Vext{^ = (l/2)m(cJ ■ r}^ a harmonic potential, 
such as in Ref.jl], and Vext{r) = \A-r\ a linear quadrupolar potential 24 1; but one can consider 



Vq/ cos (7 ■ r), or a generic 



any confining potential such as a Poschl- Teller 
power-law potentiallsi] Vext{r) = ei\x/a\'P + e2\y/h\^ + e2,\z/c\'^. We write these potentials to 



exemplify the appropriate thermodynamic variables as well as the thermodynamic limit for 
each case. 

To illustrate how the generalized variables appear and how the thermodynamic limit is to 
be taken, we shall deal here first with a classical ideal gas. Further below, we shall treat an 
interacting fluid, classical or quantum, and we shall verify the correctness of the identification 
of the variables given here. Consider, therefore, a system given by the Hamiltonian ([1]) with 
no interatomic interactions, i.e. u(|rjj|) = 0. Assume the system is in thermodynamic 
equilibrium at temperature T. The canonical partition function is. 



N / -Q 



i=l 



(2) 



where (3 = 1/kT. We assume the external confining potential to be of the form Vext = 
Vext{x/lx,y/ly, z/lz,ri) where the quantities do not necessarily have units of length and 
7] stand for other intensive parameters that we assume to remain constant throughout. 
Integration of the partition function yields 



^(T,iV,V) = ^^(C(/?)Vf , (3) 



where At = h/{2iTmkTY^'^ is de Broglie thermal wavelength, V = Ixlyh is the generalized 
volume and the function is defined by 

C{(3)V= f e-^^-*("^dV. (4) 



Helmholtz free energy is found with F = —kTlnZ and, after taking the limit oo, 
yields 



F{N, T, V) = -NkT In 



1 . (5) 



For the free energy per particle, F/N, to remain finite in the thermodynamic limit, 
oo, it must be required that the "generalized" volume diverges, i.e. V ^ oo, keeping 
constant the "density" N/V. As it will fully justified below, V is an exienszfe thermodynamic 
variable. The generalized volume certainly is proportional to the actual average volume that 
the system occupies, V ~ ({P)V. This average volume, however, is not an independent 
thermodynamic variable since it depends on the temperature. Moreover, it is not a correct 
variable since the actual volume that the system occupies is, in fact, unbounded. This can 
be seen by calculating the density profile which, as shown below, in this case is proportional 
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to 

p(f) ~ e-^^-*^'^. (6) 

Nevertheless, the thermodynamic hmit V ^ oo indeed imphes that the volume of the system 
becomes arbitrarily large. 

As it will also be justified, the generalized pressure is the intensive conjugate variable to 

V and may be calculated from 

For an ideal classical gas and for any confining potential, it is found that obeys VV = NkT, 
the ideal equation of state but in the appropriate variables. 

For the particular external potentials here considered one finds, V = V for rigid walls, 

V = 1/uJxUJyUJz for the harmonic potential, V = l/A^AyAz for the quadrupolar poten- 
tial, V = ^/'Jxlylz for the Poschl- Teller potential and V = abc for the generic power-law 
potential. Likewise, we find C{P) = 1 for rigid walls, = (27r//9m)^/^ for the har- 
monic potential, (iP) = Sn/p^ for the quadrupolar potential, CiP) = 23/3-(i/p+iA+i/g)r(l + 
l/p)r(l + l//)r(l + l/q')/ei^^e2^'e3^'' for the generic power-law potential, and ({jS) = 
4tt /^^^2 6xp[— /?Vo/ cosx]dx for the Poschl- Teller potential. As we shall find below, while 
the role of the generalized volume is completely analogous to that of the usual volume in 
homogeneous systems, the thermodynamic properties of the different confined fiuids show 
very strong variations on their temperature dependences due to the function C{l3). 

For purposes of presentation we have assumed that the fluids under study are effectively 
three-dimensional. That is, we suppose that in the three spatial dimensions the trap becomes 
macroscopic. It is clear that the theory can be adjusted to deal with (quasi) two- and one- 
dimensional systems. For this to occur, the trap must be spatially very tight in one or two 
directions, which also implies that the fluid thermal excitations in those directions are smaller 
than the trap energy levels in the same direction; thermodynamic behavior can only exist 
in the remaining directions. This is the case for the very recent studies on the Berezinskii- 
Kosterlitz-Thouless transition in quasi 2D harmonic traps in clouds of ^^Rb 33|, l3J] and ^^Na 
gases jssl- and where the generalized area and pressure can be deflned. 
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III. VIRIAL EXPANSION FOR ARBITRARY CONFINING POTENTIALS. 



With the identification of the generahzed variables and the corresponding thermodynamic 
limit in hand, we now turn to the general problem of an interacting gas, classical or quantum. 
We extend the analysis described in the texts by Mayer and Mayer |3f|, ter HaarjSTJ and 
Blatt[38]. Again, we assume the system is in thermodynamic equilibrium at temperature T 
and we analyze it in the grand canonical ensemble. We thus consider a chemical potential n 
whose value may be found by imposing a given number of particles N. The grand potential 
is given by. 



n = -kT In J2 e^^^ Tr'e"^-^^, 

N=0 



where 



if the system is classical, and 



Tr' e-^''- = ^ye' 



(f^r < fi,r2, . . .,f]sf\e 



rip, 'r2Pi ■ ■ ■ I'^NP >5 



(8) 



(9) 



(10) 



if the system is quantum. The sum is over all permutations of 1,2, and e = ±1 for 
bosons or fermions. 



To find the virial expansion, equation (IHl) is first rewritten as.|36l. 1371. l38| 

oo -| 

where the functions are given by, 

Jd^^-r Unin, ...,fn) 



(11) 



(12) 



for classical and quantum systems, respectively, and where the Ursell functions are given by 
the hierarchy: first order, 

f/i(l) = Wi(l), (13) 

second order. 



t/2(l,2) = W^2(l,2)-f/i(l)t/i(2) 



(14) 
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third order, 



[/3(1,2,3) = W^3(l,2,3)-[/i(l)[/2(2,3)-C/i(2)C/2(l,3)-t/i(3)C/2(l,2)- 

C/i(l)C/i(2)C/i(3). (15) 



and so on, and 



Wn{l,2,...,n) 



(16) 



Epe^ < fi,f2,...,fn\e 



-l3Hn 



rip,r2P, ...,rnp> 



for classical and quantum systems, respectively. 

The problem of the virial expansion reduces to find the value of each contribution in the 
thermodynamic limit, taking into account the interactions among the atoms or molecules. 
This is what we do now for a general confining external potential Vext{r)- We proceed by 
systematically calculating 7„ order by order and then generalize it to /„. We did so from Ii 
to I4. Since the calculations are quite lengthy, though straightforward, we explicitly present 
in the Appendix the case I2 only. Next we discuss the results. 

The calculation of Ji is very simple and it turns out that, in the thermodynamic limit, 
the classical and quantum cases give the same result: 



fe3 



(17) 



where the one-particle Hamiltonian Hi is given by 



In the limit, one finds. 



zm 



h = ^C(/5)v 



(18) 



(19) 



for both classical and quantum cases. 

In the Appendix we show the explicit calculation of 12- That analysis suffices to see how 
to find /„. The key is in the separation of center of mass and relative coordinates. This 
change of variables is generally, ^ = ^(ri + f2 + . . . + r„), f^^^ — ri — ^2, r^^^ = ^2 — r3, . . ., 
f^"~^^ = r„_i — f„, with their canonical conjugate momenta. 



In both the classical and quantum cases, the main assumption is the same, namely, that 
one must take the thermodynamic limit V ^ oo. This allows to make the approximation, 

where (^1,^2, . . . , r„) are to be given in terms of the variables {R, r^^\r^'^\ . . . , f^"-!)) by the 
above transformation. This approximation separates the center of mass motion from the 
relative ones. The former is always quasiclassical and its contribution to /„ is proportional to 
({nf3)V/ Xq^, while the contribution from the relative coordinates yields the virial coefficients 
bn] these can be classical or quantum, but they are the same for any external potential. That 
is, we find 

In = bn. (21) 

The validity of the above procedure, in the classical case, reduces to require that the inter- 
molecular potential must vanish for lengths r ^ a, with a the range of such a potential. 
Additionally, the interaction must be "short-range" , namely, decaying faster than 



27 



36| . In the quantum case and for high 



otherwise the virial coefficients 6„ do not exist 
temperatures, the range of the relative variables r is bounded due to presence of the potential 
u{r) and the validity has the same limitations as in the the classical case. At low tempera- 
tures the bound is set up by either the thermal de Broglie wavelength or the scattering length 
a 391]. If the gas behaves as an ideal one, the relevant length is de Broglie wavelength. In any 
case, as long as the relative coordinates remain bounded by a finite quantity, however large, 
one can take the limit of very large volumes, V ^ 00, and implement the thermodynamic 
limit just as in the classical case. 

Summarizing, we find that in the thermodynamic limit the grand potential can be written 
in general as, 

^=-kTVY.-r-^UT) (22) 

n=l "'■ 

where quantum or classical virial coefficients should be used; in both cases 61 = 1. Expression 
fl22|) for the grand potential is one of the main results of this article. It is the virial expansion 
for arbitrary confining potentials. The number of particles and the entropy S can be 
calculated from (minus) the partial derivatives of Q with respect to fi and T respectively. 
Q, N and S are found to be homogeneous first order functions of V, and this implies that V 
must be an extensive variable and justifies the thermodynamic limit as used above. From 
the expression fl = —W, the generalized pressure is read off (!22l) . 



Note that the most important difference of the grand potential between a given arbitrary 
external potential and the homogeneous case is the function C{l3) rather than the generalized 
volume V. The latter enters in the same way for any potential, including the rigid-walls case; 
that is, it gives rise to the intensive quantities formed between the extensive variables A^, 
S, E etc. and V, and that remain finite in the thermodynamic limit, i.e. N/V, S/V, E/V, 
etc. However, as it is well known from calculations in external potentials, the temperature 
dependence of the thermodynamic variables is very different and unique for each external 
potential. This difference is contained in the function C{f3). This will be exemplified in the 
next section. 

To illustrate the use of Eq. fl22l) . we apply it to an ideal quantum gas. From the analysis 
in the Appendix and their corresponding value for third and fourth orders, one finds that 
the virial coefficients are given by 

&r = e"-^:iAr-^). (23) 



n 



The grand potential for an ideal quantum gas can, thus, be written as 

-/5^=3TEe"^^C(n/5)-^. (24) 

'^T n=l ^ 

This formula can be directly compared with the corresponding ones for, say, the rigid 
walls potential V = V and C{n(3) = 1, or the harmonic potential V = l/u^ and 
C{nl3) = (2TxkT /nmY^'^ . The "textbook" formulae for these potentials are, [4^ 

y 1 x^/^rfx 



for rigid walls, and 

for a 3D isotropic harmonic potential. Expansion of the integrals of these last two equations 
in powers of e^^ yield the virial expansion, Eq.( 



IV. PHYSICAL MEANING OF THE GENERALIZED PRESSURE. EQUATION 
OF STATE AND HEAT CAPACITY. 

In addition to the previous discussions, we want to emphasize several points that we 
consider to be of importance for both the understanding of the present results and their 
practical applications. 
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First, we recall a result previously set forth by one of the authors |15l. Il6|, [iTj. This is the 
fact that since V is an extensive generalized volume, there exists an intensive generalized 
pressure V, conjugate to V, given by 

and, therefore, Q = —VV. A relevant consequence is that Euler relationship is = TS — 
VV + for any external potential in the limit here discussed. Thus, V being a bona- 
fide thermodynamic variable allows us to calculate, for instance, the Gibbs potential G = 
G{T, N,V) or the entalphy H = H{S, N,V) by appropriate Legendre transforms. The 
identification of the generalized pressure also yields the correct equation of state of the fluid; 
this is the relationship V = V(y/N,T). 

Second, the identification of V is not only a formal one but it has a clear physical in- 
terpretation as the variable responsible for mechanical equilibrium. It must be recognized 
that the usual hydrostatic pressure is not the appropriate thermodynamic variable in an 
inhomogeneous fluid since it is a local variable. In general, however, mechanical equilibrium 
of an inhomogeneous fluid is given by Pascal law, 

V-P(f) = -p(f)VV;.t(r1, (28) 

where P(r) is the pressure tensor of the fluid and p(r) is the local particle density, namely, 
the density profile. As we shall see below, barring phase-separated states within the confined 
fluid, one expects the pressure tensor to be a local quantity, P{r) = p{r)l, where 1 is the 
unit tensor and p{f) the local hydrostatic pressure. Integration of the virial of the external 
force, namely, of — r ■ {p{r)'Wext{'>^) , yields 

J Tip (fr = J (fr pif) r ■ Wext{r) 

= l^r,-y^ye.t{n)Y (29) 

where in the second line we have identified the average of the virial of the external force. 
This average can be performed in any ensemble. Let us choose the canonical ensemble. The 
partition function is 

Z = Tr e'^^^ (30) 

with H]^ the Hamiltonian of the interacting particles confined by Vexti'f^-, as given by 
Eq.([T]). Recalling that Helmoltz free energy is F = — kT In Z, one finds F = F[N,T,V), 
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and by a simple calculation, 

V 



'dF 



T,N 

= ^l^n■V^VUr^)Y (31) 

This result has a two- fold relevance. Comparing with Eg. (1291) . one finds the physical signif- 
icance of the generalized pressure, 

PV =- I TtP (fr. (32) 
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It is simply an intensive global property of the pressure tensor, the quantity that bears the 
information of mechanical equilibrium in the body. An alternative form is clearly, 

= i 1 p{r) f ■ V\/ext(r) (33) 

which may be seen as the generalization of the celebrated theorem found by Clausius in the 
middle of the XIX century for homogenous systems: The (average) virial of the external 
force equals —?)W. It follows as well from the present discussion that the reversible work 
performed on or by the system after a change of the external potential is given by, 

dW = -V dV. (34) 

Further, in addition to the interpretation of the generalized pressure, we point out that 
Eq. (!33!) is a remarkable formula since it gives rise, with the sole knowledge of the density 
profile p(r), to a direct route for the calculation of the generalized pressure and, hence, to 
the equation of state of the fiuid, V = V{V/N,T). We recall that the density profile may be 
known from exact or approximate calculations, or directly from experiments. It is somewhat 
puzzling to realize that expression (|33|) yields only an identity for the rigid-wall case and does 



not give a calculational tool for the hydrostatic pressure. The latter needs the knowledge, at 
least for pairwise interatomic interactions, of the two-body density correlation function [4 li]. 
Here, we find that for inhomogeneous systems knowledge of the one-body density suffices. 



This was the procedure used in Ref.[l5l] to obtain the equation of state of a degenerate 
interacting confined Bose gas within the Hartree-Fock approximation. 

For completeness of our presentation, we write down the first few terms of the so-called 
virial expansion of the equation of state V = V{N/V,T). This can be done by finding 

12 



= N{fi,T,V) from Eq. (l22|) and inverting it term by term to yield fi = fi{N/V,T), then, 
substituting the result into V = V{fi,T): 



N 



N 
V' 



kT 



2 ^2(^) 21 



3 cm 



+ ... 



(35) 



Once again, we remark that the functions Cil^) make all the difference. Since in some 
instances one can refer the calculation of the virial coefficients to a diagramatic expansion 361] , 
one finds that the diagrams sum up differently for different potentials. We also recall that 
this type of virial expansion was used in Ref. 2J| to fit experimental data from a gas of 
Sodium atoms in a quadrupolar potential. 

As a further instance of the relevance of these variables, we mention here that there has 
been a lot of interest recently on the behavior of fermi gases (^°K and ^Li) near the unitarity 



limit where the scattering length diverges, because it appears that thermodynamics 



universal, i.e. independent of the interatomic interactions; see e.g. Refs 
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becomes 



Id, 
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231]. In this region, the gas behaves as an ideal one in the sense that it obeys the ideal 



virial relationship E = 2N < V^xt > for a harmonic trap, with E the internal energy. The 
relevance to the present work is that < Vext > for a harmonic trap is 3/2'PV, whether 
the gas behaves as ideal or n ot, see Eq. (!33|) and Ref.[l5|. That is, the quantity that has 



been measured in Refs 



.0, 0,y, 



231 ] is precisely the generalized pressure for the harmonic 



trap. And indeed, if the universality hipothesis is correct 211], using the virial theorem for 
ideal gases for arbitrary potentials, yields the following relation that should be obeyed in 
the unitarity region, 

N q 

(36) 



N o 

E=<Y.Vexm>+-VV. 

1=1 ^ 



Away from the unitarity limit, this equation is no longer valid but the measurement of V, and 
of the heat capacity as we describe below, can be performed to obtain the thermodynamics 
of those states. 

The equation of state does not suffice to know all the thermodynamic properties of a 
one-component fluid. One needs a second relationship such as the heat capacity at constant 
generalized volume Cy = C\;{N,V,T). We propose now a very simple measurement to 
find out such a heat capacity in the current ultracold fluids. This consists of an adiabatic 
compression or expansion: we note that the ultracold trapped gases are actually isolated 
and confined by magnetic or optical traps [3, y, 0, ^, sl, Q, 



101, 



ll| and, therefore, a 
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slow change of the confining potential, namely, of the generalized volume, will give rise to 
an increase or decrease of the temperature, depending on whether the fluid is compressed or 
expanded. This procedure is adiabatic. Since the generalized volume and temperature are 
measured in the current experiments, the quantity ((9V/9T)jv,s can thus also be measured. 
Then, from a thermodynamic relationship one can evaluate the corresponding heat capacity, 

'-'(S)„(S)„ 

where previous knowledge of the equation of state is needed for the calculation of the second 
factor on the right hand side. Note that the measurements of the equation of state and of 
the heat capacity correspond to two different sets of experiments. 

V. A COMMENT ON THE "LOCAL DENSITY APPROXIMATION". 

The exactness, or validity, of the "local density approximation" follows right away from 
the corresponding expressions for Q, N and S, as given by Eq. fl22|) . Consider the rigid- wall 
external potential, V = V and C{n(3) = 1. We can thus deflne the grand potential per 
unit volume u){fi,T) = ^l/V, the number of particles per unit volume (particle density) 
p{fi,T) = N/V and the entropy per unit volume s{fi,T) = S/V. We now consider the 
same system but conflned by an external potential Vexti'f^- We can flnd its thermodynamic 
properties by implementing the "local density approximation": take uj, p and s of the 
homogenous case and make those functions per unit volume to be their "local" densities 
co'(r), p{r) and s(r*), in the presence of the given external potential, by replacing the chemical 
potential p by the "local" chemical potential piocaiif) = ~ Vextif)- It turns out that 
integration of co'(r), p{r) and s(r) over all space yield the exact expansions for f2, and S, 
in the presence of Vext(^), as given by Eq.( |22l) and its derivatives. That is, one flnds that 
LDA procedure gives rise to exact results. We recall that the validity of LDA for classical and 



quantum systems in this limit was rigorously proved in Refs. [27|, |28|] and 29|] , respectively. 

The above description does show that in the thermodynamic limit the system is locally 
homogenous and that "locally" actually means in length scales large compared with those 
of interatomic interactions. It is in this latter connection that LDA is largely used without 
the need of further justiflcation. There is a warning, however, that must be raised when 
using LDA. It may appear that if one is able to flnd any thermodynamic variable q for a 
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homogenous system and express it in terms of the chemical potential and temperature 
T, namely q = q{fi,T), its local counterpart when in the presence of an external potential 
Vexti^ is simply g(r) = q{fiiocai{^,T). This, in general, is incorrect; its is strictly justified 
for ^l/V, N/V and S/V only . It is incorrect, for instance, for the internal energy and other 
free energies, except Q, as well as for other functions such as heat capacities. This statement 
can be verified by using expression (122!) for the Grand Potential of a confined fluid. Thus, 
the fact that the system is locally homogeneous does not imply that the local states of 
the confined fluid are thermodynamic states of the corresponding homogeneous system, i.e. 
qiocai{r) ^ q{lJLiocai{r),T) in general. 



VI. FINAL REMARKS 



As a final comment, we emphasize that the use of the generalized thermodynamic vari- 
ables here presented should lead to correct results by simply following the rules of thermo- 
dynamics, without necessarily resorting to a local picture. It is clear to us that for this to 
be useful one needs, first, to change the usual intuition on volume and hydrostatic pressure 
to their generalized counterparts and, second, to provide examples where these variables 
lead to new insights. There are already few studies where the latter has been put to use 
in approximated calculations of equations of state 15|, ll7|] and there already incipient exper- 



iments that have made use of them 



24 



25 



261]. Our contention is that the measurement 



of the equation of state and the heat capacity in terms of the generalized thermodynamic 
variables should offer a complementary and useful tool for the analysis of ultracold trapped 
fluids. 

Although much has been learned with use of LDA, we have indicated that care must 
be taken when using it. This should be more notorious when dealing with phase-separated 
fluids where it is not clear if LDA suffices for their description since the interfacial widths 
of the phase boundaries are expected to be of the order of the range of the intermolecular 
lis situation appears to be the case for the states found recently in trapped 



41|]. T 



mteractions 
gases of ®Li atoms 



where it is found that the confined fiuid phase-separates into a 



superfiuid and a normal paramagnetic gas, showing clearly an interfacial phase boundary. 
There are differing theoretical studies of whether LDA should be enough or if "surface 



tension" terms should be included 



18 



,[l9 



In general, for such inhomogeneous states. 



15 



one should not expect a local picture to be valid across the interface; the thermodynamic 
potentials are indeed expected to be non-local on the density profiles. 
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APPENDIX: DERIVATION OF h 

In this Appendix we provide the derivation of the second order contribution to the virial 
expansion. 

Classical case. From Section III we find that I2 in the classical case is given by 

1 (I rf^pi j d'n e-^^^) (/ d'p, I d'r, e-'^^^) (A.l) 

with H2 given by eq.([T]) for N = 2. We make the change of variables to center of mass and 
relative coordinates, 

R={ri + f2)/2 and r = ri — r2, (^-2) 
with their canonical momenta P and p. Thus, 

^Jd'Rj dVe-^(^-*(^+5)+^-(^-i^» (A.3) 
where the center of mass Hamiltonian is 



p2 



and the 2-particle relative-coordinates one is, 

2(4) '^-^ 

The integrals over the momenta yield 1/Aj^. Notice, however, that the integrals over R 
depend strongly on r, even for particles confined in a rigid vessel of volume V . In such a 

16 



case, the boundary terms couple the integrals. We rewrite I2 as follows, 

1^ = 1^ j d^Rj rf3^g-;3(V...(i?+£)+K..(i?-£)) (A.6) 

where we have introduced the Mayer function, /(r) = e~^'^^^^ — 1. Here is the important 
step: the function /(r) vanishes for distances r longer than the range of the intermolecular 
potential a, i.e. the value of r is bounded, r < a. Hence, in the thermodynamic limit, in 
which the external potential becomes "shallower" and "shallower" and the volume of the 
system larger and larger, we can set, 

— * — * 

VeAR+l) + ^ 2Ve.tiR). (A.7) 

This is the thermodynamic limit. Therefore, 

l2 = ^ [ d'Re-'^''-^^^^ I d'rfir) = ^Cmh, (A.8) 

where we have identified the classical second virial coefficient 62- Thus, the well-known 
bound is that the intermolecular potential must vanish for lengths r ^ a. Clearly, it must 
be "short-range" interaction (faster than 1/r^), otherwise 62 does not exist. 
Quantum case. Again, from the expressions in Section HI, one finds, 

I2 = j d^'^i j ^^^"^ ri,r2|e"^^^|fi,r2 > +e < fi, r2|e"^^^|f2, fi - 

j Sri < ri|e"^^i|fi > j d?r2 < f2|e"'^^^ |f2 > (A.9) 

where e = ±1 for bosons or fermions. The derivation follows essentially the same steps as 
in the classical case. First, one performs the same change of variables as in Eq. (]A.2l) to 
rewrite J2. Then, the thermodynamic limit may be taken by separating the center of mass 
motion from the relative one. It can be realized that for high temperatures the range of the 
variable r*is bounded due to presence of the potential u{r). At low temperatures the bound 
is set up by either the thermal de Broglie wavelength or the scattering length a. Therefore, 
as long as r remains bounded by a finite quantity, however large, one can take the limit of 
very large volumes, V — > 00, and implement the thermodynamic limit just as in the classical 
case. One finds, 

h = ^cmK (A.io) 
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formally identical with its classical counterpart, formula (lA.Sp . but with the quantum second 
virial coefficient: 



62 = 2^/2 A| J d^r (< r\e 



'I3H2 



\f> + 



r > 



< r\e 



-/3pV2{m/2) 



r > 



(A.ll) 



where 0,2 is given by Eg. (]A.5[ 



motion is always quasiclassical 



. As a rule, in the thermodynamic limit the center of mass 
421 ] ■ The expression for the quantum second virial coefficient 
above can be seen to be the correct one by comparing, for instance, with the expression 
given in Ref. 



42| . For slow collisions, the relevant ones for ultracold gases, 62 depends on the 
scattering length and this may become quite large near a Feschbach or potential resonance. 
The formulae here derived may then be not applicable very near such a point, called the 



unitarity limit, but as it has been shown 



2l[ | this may be expected since in such a limit the 



system behaves as if near a critical point. We add that the description of the scattering al 
low energies near resonances is valid for interatomic potentials u{r) that decay at least as 
l/r3[39|. 

To end this part, we find illustrative to calculate 62 for an ideal quantum gas, i.e. for 
u{r) =0. One finds the so-called "exchange" contribution to the second virial coefficient: 



l(0) 



A3 

23/2 ^T- 



(A.12) 
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